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Computations of the strong field generation of gravitational waves by black hole processes produce waveforms that 
are dominated by quasinormal (QN) ringing, a damped oscillation characteristic of the black hole. We describe here 
the mathematical problem of quantifying the QN content of the waveforms generated. This is done in several steps: 
(i) We develop the mathematics of QN systems that are complete (in a sense to be defined) and show that there 
is a quantity, the "excitation coefficient," that appears to have the properties needed to quantify QN content, (ii) 
We show that incomplete systems can (at least sometimes) be converted to physically equivalent complete systems. 
Most notably, we give a rigorous proof of completeness for a specific modified model problem, (iii) We evaluate 
on : the excitation coefficient for the model problem, and demonstrate that the excitation coefficient is of limited utility. 
We finish by discussing the general question of quantification of QN excitations, and offer a few speculations about 
unavoidable differences between normal mode and QN systems. 
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I. INTRODUCTION AND OUTLINE 

Essentially all computations of the generation of gravitational waves by strong field black hole processes produce a 
gravitational wave with the shape of a damped sinusoid jlj . The oscillation period and damping time depend only on 
the parameters of the black hole, and not on the manner of excitation. The meaning of the complex frequency of this 
damped oscillation is now well understood. A single frequency perturbation outside the hole can satisfy the natural 
radiative boundary conditions (radiation into the black hole and outward to infinity) only if the frequency is one of 
the discrete set of frequencies, called quasinormal (hereafter QN) frequencies. The least damped of these complex 
frequencies is what dominates the appearance of computed waveforms. 

QN excitations are relevant, in principle, to most or all systems with radiative boundary conditions. Stellar models 
for example have short periods for nonradial oscillations driven by fluid pressures, and long damping times of these 
fluid oscillations due to the weak emission of gravitational waves. The motions of the stellar fluid can be studied with 
radiation damping omitted (e.g., with the use of Newtonian gravitation theory, or Post-Newtonian theory) and the 
weak radiation can be added, after the fact. When the radiative coupling is "turned off" the problem of the oscillation 
of a perfect fluid stellar model can be analyzed in normal modes and one can find the radiated energy coming 
from each separate oscillation frequency, and can decompose the radiative power into that fraction assigned to each 
frequency. 

The situation is dramatically different for black holes, which have only a single time scale. (For a nonrotating hole 
this is 2GM/c 3 where G is the universal gravitational constant, c is the speed of light, and M is the mass of the 
hole.) The period and damping time are therefore of the same order and there is no meaningful way of turning off 
the damping for black hole oscillations; there is no underlying normal mode system. This suggests that there may 
be no clear way of specifying "how much QN ringing" of some particular black hole QN frequency is contained in an 
emitted waveform. This suggestion is made plausible by the mathematical origins of normal modes and QN modes. 
The properties and usefulness of normal modes are closely related to the fact that they are eigensolutions to a self 
adjoint problem. QN modes, on the other hand, are eigensolutions of a problem that is not self adjoint. But the 
dominance of QN frequencies in computed waveforms is so robust that it seems that the strength of QN ringing must 
be quantifiable, or at least that mathematical sense must be made of the question. 

In this paper we try to make mathematical sense of quantification. In attempting this we draw upon parallels with 
normal modes systems. By the "excitation" of a mode we mean, in parallel to excitation in normal modes systems, an 
index of the contribution that each mode makes to the overall waveform and to the energy. To develop a description 
of QN excitation we start with a viewpoint that a meaningful and rigorous quantification is very implausible unless 
the QN system is, in some sense, complete. We then follow a three step process. First, in Sec. II, we define and posit 
the existence of QN systems that are complete (in a sense to be defined) . We then point out difficulties in quantifying 
excitation in a complete QN system. We construct a particular measure, the "excitation coefficient" that overcomes 
these difficulties, and is closely related to the description of the excitation of normal modes. 

Our next step is to prove the existence of complete QN systems and relate the mathematics of black hole processes 
to complete systems. This step, carried out in Sec. Ill requires a rather lengthy discussion of "induced completeness." 
Though this discussion is not directly related to the problem of QN excitation, it is a necessary step (and is interesting 
in its own right). The discussion in Sec. Ill shows that completeness can be induced. That is, an incomplete QN 
system can be changed with a modification that satisfies two criteria: (i) The effect of the modification can be made 
arbitrarily weak. More specifically, the modification can be made small enough so that the waveform that evolves 
from any initial conditions is arbitrarily close to the waveform evolving with no modification, (ii) No matter how 
weak the modification is, the modified QN system is complete. Our demonstration in Sec. Ill does not consists of 
a general theory for such modifications; a conjecture about the general conditions has been given by Young et al. 
H (though their definition of completeness is somewhat different from ours). Here we will sacrifice generality and 
direct astrophysical relevance for specificity and rigor. We present the details of a specific model. We will start with 
a model, the 'TDP' with only a single conjugate pair of QN frequencies, and modify it to the "spiked TDP," a model 
with an infinite QN spectrum. The Appendix gives a rigorous proof of completeness of the spiked TDP (i.e., that 
under specified circumstances the outgoing waveform is a convergent sum of components at quasinormal frequencies). 
Numerical results are shown in Sec. Ill to demonstrate the negligible effect of the modification, and to demonstrate 
the pattern of convergence of sums of single frequency excitations. 

Having established that completeness can be induced (at least in one model problem), we return, in Sec. IV, to 
the question of measuring the excitation of QN modes and in particular to the excitation coefficient, introduced in 
Sec. II. We demonstrate, with a few examples, that this formal measure of excitation does not generally give a useful 
quantification. The failure of this measure is discussed along with a broader discussion of differences between QN and 
normal mode systems, and conjectures about mathematical properties of QN systems. 
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II. COMPLETE QN SYSTEMS 

A. Definition of QN frequencies 

For definiteness we will limit considerations to solutions of the equation 

9 2 * d 2 * 



dx 2 dt 2 



V{x)^ = . (1) 



Such an equation describes the dynamics of many mechanical systems, and the evolution of multipole perturbations 
(scalar, electromagnetic, or gravitational) of spherically symmetric (Schwarzschild) black holes |§-|7]. Perturbations 
of rotating (Kerr) holes 0, on the other hand, cannot be reduced to radial-time equations. QN oscil 
frequency solutions of the form r) = ip(x) exp (itot) and hence are solutions of the equation 



ations are single 



+ [u 2 V(x)] * = . (2) 

We assume that the domain of x includes x — oo, and that the nature of the potential V(x) is such that a "radiative 
boundary condition" can be defined at x — > oo. A clear example is a potential with support of V(x) only for x less 
than some x max . In this case the boundary condition is that ij)(x) oc exp (— ilux) for x > x ma x- 

For potentials that do not vanish, but fall off sufficiently fast as x — > oo, the more general radiative boundary 
condition for real frequencies will be that for large x, the solution for ip(x) have the form exp (— iu>x)F(u>, x), with 
F -^constant, as x — > oo. This condition is not quite sufficient if u> has a positive imaginary part; see [|8| for a 
complete discussion. In short, the solution satisfying a radiative boundary condition for complex u) can be regarded 
as an analytical continuation of a solution satisfying a radiative boundary condition for real to. 

The boundary condition at the other end of of the x domain may be a standard Sturm-Liouville boundary condition 
(e.g., 'J = at x = 0) or may be a radiative boundary condition at x — > — oo. For spherically symmetric black holes, 
the range of x extends from — oo to oo, with — oo representing the black hole horizon The potentials fall off 

exponentially in x as x — > — oo and as const/x 2 , as x — > oo. Radiative boundary conditions are imposed both at 
x — > — oo and x — > oo, corresponding to radiation moving inward through the black hole horizon, and outward towards 
spatial infinity. 

QN frequencies, are the eigenvalues to = loqn to the problem defined by (^|) for radiative boundary conditions of 
the type just discussed. Due to the boundary conditions, this problem is generally not of the Sturm-Liouville type 
and the usual features of eigenvalues of a Sturm-Liouville problem are absent. In particular the QN frequencies loqn 
are generally not real. A positive imaginary part indicates an exponential decrease with time. A negative imaginary 
part would indicate an instability; no frequencies with negative imaginary parts have been found for black hole QN 
systems. 

It is clear that QN frequencies must occur in conjugate pairs. If loqn is a solution to the eigenproblem corresponding 
to tpQN, then —luq N is also a solution corresponding to iPqn- We will use a tilde (~) to denote the conjugate relationship 
of QN frequencies. Thus the conjugate to QN frequency co^ is u>j, that is loj = —loj. 



B. Definition of completeness 

We will choose to give a rather specific meaning to a complete system of QN modes. The Cauchy data for (|l|) 
consists of ipo and ipo, the initial value of $ and of its time derivative, 

M*) = *(*.s)lt=o Mx) = d^(t,x)/dt\ t=0 (3) 

We consider an interval X2 < x < x\ in the domain of x and we consider Cauchy data at t = for ([l]) which has 
support only in this interval. We then consider the solutions $ to ([!]) for such data, and we focus attention on the 
value of this solution at x b s , a particular value of x satisfying a; bs > x\. This corresponds to the physical situation 
of an observer at a; b s detecting radiation resulting from an initial disturbance located at some distance from her. 
The "observed waveform" that we focus on is then 



f{t) = ^!(x ohs ,t) 



(4) 
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FIG. 1. Propagation of initial data to the observation location x b s defining the waveform f(t). 

As shown in Fig. ^, there will in general be a minimum value t m - ln of i, such that the point x Q \> s ,t is influenced 
by the evolved Cauchy data. That is, for t < t min , the area between the past directed characteristics from t,x ^ s 
intersects the t — hypersurface outside the support of the Cauchy data. We are interested in fit) only for t > t m i n . 
The physical interpretation of this is that we are considering only the waveform generated by the Cauchy data. 

We take a complete QN system to be one which satisfies the following criteria: 

• The solutions to the QN eigenvalue problem form a discrete spectrum and can be arranged in order of increasing 

• We consider only Cauchy data that 

— has support only within a compact region [2:2, ati]. 

— belongs to a specific continuity class C p , where p depends on the nature of the problem. 

— results in a wave form which is square integrable from t — t m i n to oo. 

• For such Cauchy data, the waveform fit) that evolves from any such allowed Cauchy data can be written as 



/(*) 



(5) 



Here a™ is the n th coefficient in the sum over QN modes. Since f(t) is a function of x bs, the a™ coefficients are 
also functions of Xobs, but we shall not explicitly exhibit this dependence. The summation in ([5]) is in order of 
increasing |K(w„)|, and the convergence is uniform for t > i m ; n . 

It is important to note that our view of completeness is rather different from other possible meanings of the term. In 
particular, our choice of the meaning of completeness has nothing directly to do with the x-dependence of the single 
frequency solutions and with the question of whether these solutions can be used to span acceptable Cauchy data. 
Our meaning of completeness, then, is rather different from that of Young et al. jj]. It also disagrees with the concept 
of completeness used by Husain and Price and by Beyer and Beyer and Schmidt Q 



C. Function space and inner product 

In accordance with our definition of completeness, our vector space is the space of all functions f(t),t > £ m j n that 
can evolve from acceptable Cauchy data. Our class of acceptable Cauchy data will always be chosen so that fit) is 
square integrable from t = t mlu to oo. On this space of functions we define an inner product to be: 

/>oo 

f-9 = J t r(t)g(t)dt. (6) 

We could of course include a weight function W(t — t m { n ) in the integral defining the inner product, but the time 
translational symmetry of the background suggests that W should be constant. The choice in Eq. (||), furthermore, 
means that / • / is the time integral of the square of the wave function, a measure closely related to the energy content 
of a wave. (For black hole processes, the connection with gravitational wave power will be made explicit presently.) 
Our assumption of completeness above means that the functions exp (iwQjvi), while not elements of our function space 
themselves, span this function space in the sense of (Q); we will therefore consider them a basis. For a function f(t) 
in our space we can use the inner product to compute another set of coefficients a„ by 
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a„ = (e^*) • /(t) = jf e-^f(t)dt. (7) 
The following relations for the coefficients of conjugate modes are straightforward to verify: 

H = («*)* « s = (« fc r (8) 

Since the convergence is uniform by hypothesis, we can integrate term by term in the sum-of-modes expression for 
the norm of / to find: 



J t °° (X>" e ^*) /(*)* = X> n )* 



\f{t)?dt = / j yy <'-■■' \ /u),n = >>"> <•„ . (<)] 

The final sum in (||) is real, as it must be, since for any k, the sum (a k )*ak+ (a k )*aj c is real. 

In most physical problems the radiated power is the square of the time derivative of the waveform. If at x b s 
our wave function evolving from the initial data is f(t) = ip(x i, s ,t). If /(i m in) vanishes (i.e., if the waveform starts 
continuously) then this type of energy can be evaluated as 

\f(t)\ 2 dt = J2«f(a n ran. (10) 

» n 

As in (||) the reality of the sum is guaranteed by the relations of conjugate coefficients in (@). 

Since we have an inner product, we have an equivalence between vectors and dual vectors in our function space and 
we can define a set of covariant basis functions (j> m (t) by the property (j> m (t) ■ e M ™* = S nm . If follows from Eq. (Q) that 
the a n are the expansion coefficients for f(t) with respect to the covarariant basis functions <f) n . We shall henceforth 
refer to a n and a n , respectively, as the contravariant and covariant coefficients of fit). The components of the metric, 
in this function space, with respect to the QN basis, are 

y^-Odt = c nk . (ii) 

It should be noted that G is a hermitian matrix, but it is not diagonal, i.e., the QN oscillations are not orthogonal. 
The metric coefficients can be used, in principle, to relate a n and a n . The expression for f(t) in Eq. (^) can be 
substituted in Eq. (Q). Since the convergence in Eq. (|^) is uniform, we can integrate term by term and get 



5>. 



k a k . (12) 



Since G n k is not diagonal, the covariant basis vectors, i.e., the basis vectors dual to e 1 ""*, are mixtures of the e luJnt 
functions (usually involving all of them) . An indication of the unfamiliar problems this produces can be seen in the 
following rough argument. Let us suppose that we have a waveform that in some sense is "almost pure" (say) seventh 
mode. That is, suppose that f(t) » a 7 e 1 " 7 * + a 7 e lu "' t . (We are supposing that this relationship is only approximate 
since it will in general be impossible to excite a truly pure single frequency mode with smooth, compact initial data.) 
This waveform, for which (almost) the only contravariant coefficients are a 7 and a 7 , will have contravariant coefficients 
a n for all n. A waveform that is a pure (or almost pure) single mode excitation in one sense is therefore not a single 
mode excitation in another. This presages some of the problems in quantifying the excitation of a mode, and we will 
return to this point at the end of Sec. IV. 

It is possible, of course, to use Gram-Schmidt orthogonalization to find basis function which are orthogonal according 
to the inner product of (^|). The resulting basis functions will not (except for one of them) correspond to single 
frequency excitations, and do not seem to be of interest. 



D. Intuitive insights; the excitation coefficient 



Some rough considerations of model problems suggest the intuitive basis for some of the mathematical difficulties 
to appear below, and point to a possible approach to quantifying excitation. The most obvious difficulty is the 'time 



shift' problem. We imagine a configuration like that pictured in Fig. ( 
two sets of Cauchy data (a) and (b). The two Cauchy data sets are loca 



3): A potential with compact support and 
ized and are identical except that set (b) is 
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shifted to the left, to smaller x, by some finite displacement Ax. The support of neither Cauchy data set overlaps 
the support of the potential. In this case it is clear that the waveform generated from the two cases will be identical 
except that the waveform from (b) will be shifted to later times, relative to that for (a), by an amount At = Ax. Each 
contravariant coefficient will therefore be larger in the (b) waveform than in the (a) waveform. The the conjugate 
pair of contravariant coefficients {a 7 , a 7 } will, for example, be larger by exp [(Qu^At] for (0) than for (a), though 
the excitation is physically identical. The analogous difficulty does not arise for normal modes; since they are not 
damped, the time delay only causes a phase shift. The trend is opposite for the covariant coefficients {07, a 7 }. As 
defined in (0) {07, ay} will be smaller for (6), since exp (—iui^t) (or exp (—ioj~t)) are smaller by exp [(—3^7) At] at the 
later times during which the (b) waveform has support. Therefore, neither the coefficients a k nor a k alone can provide 
a useful measure for the excitation of a QN mode. 

Any useful measure of excitation must give the same result for the two waveforms in Fig. (Q). We take advantage 
of the opposite tendencies of the contravariant and covariant coefficients under time shift to define a quantity that is 
the same for both waveforms, the "excitation coefficient" A k for the k th QN mode: 



A k = (a k )*a k + {a~ k )*a k 



(13) 



We conjecture that these excitation coefficients, or quantities constructed from them (sums of excitation coefficients, 
functions of excitation coefficients, etc.) or quantities very closely related (see the energy excitation coefficient, below) 
are the only relevant mathematical objects in the vector space that are unaffected by the time shifts, and have the 
additional properties that we outline in the following. 

In addition to the insensitivity of the excitation coefficient to time shifts of waveforms, the excitation coefficient 
has another important and relevant property (cf. Eq. (j^)): 



\f(t)\ 2 dt = J2A k 



(14) 



where the sum is over conjugate pairs. The excitation coefficients of the complete set of QN modes sums to the norm 
of the waveform. We can define a quantity closely related to A k : 



E k ee {uj* k ) 2 {a k ya k + K) 2 K)V 



(15) 



We shall refer to the E k as the energy excitation coefficient. According to (|To|), the sum over conjugate pairs of these 
coefficients gives the norm of /. The summation properties of the A k and the E k are important in clarifying what we 
mean by the "excitation" we are attempting to quantify. These coefficients appear to tell us the contribution made 
by each mode to a measure of the waveform. In the case of black hole perturbations it turns out that it is possible 
to make an even more direct connection. If f(t) is a solution of the Zerilli equation |(| for even parity perturbations, 
then the integral on the left of (jl 7 ^) is proportional to the radiated energy and E k has the appearance of the energy in 
the k th mode. If f(t) is a solution of the Regge- Wheeler equation || for odd parity perturbations, then the integral 
in Eq. (m is the energy and A k has the appearance of the energy in the k th more. 



potential 

V(x) 




Cauchy data (a) 
Cauchy data(b) 




FIG. 2. Initial data shifted in location produces identical waveforms shifted in time. 



The possibility of quantification with the excitation coefficient (or energy excitation coefficient) will be a central 
focus, of Sec. IV, but before we begin specific computations, there are a few more possibly useful insights that can 
be found from intuitive considerations. For one thing, it is interesting that the two sets of coefficients can be related 
to two different aspects of an emitted wave. The contravariant coefficients, telling us how much of a certain mode 
must be added in order to get the waveform, can be considered a "theoretical" coefficient. For a given waveform, 



Quantifying QN excitations 



7 



the projection operation on the waveform defined by (0) can be considered to give the "experimental" excitation 
coefficient. 

Simple considerations can also produce insights into the nature of the metric matrix G n k defined in ( |ll| ) . Consider 
the signal produced by some Cauchy data, and let the vector a denote the contravariant coefficients of the resulting 
wave form with respect to the QN mode basis. Let the vector of covariant coefficients be given by a, so that a = Ga. 
Now consider the same Cauchy data, shifted to the left by Ax. The new contravariant and covariant coefficient vectors 
will be given by: a k = exp (— iw^Ax) a k and a k = expiwjJAx) a k . Taking the standard linear algebra norm,, we have 
||§|| = a||o|| and ||o|| = /3||a||, where a > 1 and < 1 . 

The exact magnitude of a depends on the distribution of the contravariant coefficients and the shift Ax. In fact, 
we can always make the shift large enough so that a 3> 1. In the same way, we can ensure that (3 -C 1. Assuming 
that the metric matrix G has a minimum eigenvalue e m i n and a maximum eigenvalue e max , we then find: 



< 



L 

a 



< 



< e n 



(16) 



Since G is an infinite dimensional matrix there need not be a maxi mum or a minimum eigenvalue. For the model 
problem introduced in Sec. [II, we will give numerical evidence, in Sec. IV B, that the ratio of maximum and minimum 
eigenvalues diverges. The G matrix therefore is in some sense singular. 

We now make the additional assumption that we have constructed Cauchy data such that the wave form consists 
almost exclusively of ringing in a single QN mode. (In the specific example discussed in the following sections,, this 
will be possible for the fundamental mode using an arbitrarily short "burst" of initial data.) Let us modify the 
example of Fig. 0, by reducing the size of Cauchy data (b) by a factor exp (3wi)Ax, so that once the response to the 
(b) data starts it is identical to that of the (a) response at the same time. This situation is shown in Fig. ^ in which 
the response f(t) is the same in both plots for t > t s = At. 



potential V(x) 
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FIG. 3. Oscillations produced by shifting and scaling initial data. 



Now let us suppose that the waveforms are dominated by ringing of the fundamental mode, and that f(t) consists 
only of ringing of the fundamental mode. The difference between the two responses is a sum of QN modes which has 
the appearance of truncated fundamental mode ringing, but which has a 1 — a 1 = in its mode sum. The time t s 
can in principle be made arbitrarily long, so the truncation in the difference curve can be made arbitrarily late. We 
have then a sum of modes that looks arbitrarily close to fundamental mode ringing, but containing no fundamental 
mode. This suggests that the fundamental mode can "almost" be built as a superposition of the modes other than 
the fundamental mode, and that the modes are "almost" linearly dependent. We will discuss this further, based on 



numerical results, in Sec. [VC 



Again, this suggests that the infinite metric matrix G n k is in some sense "almost" singular. A physical insight can 
be associated with this nature of the metric: Since G n k is "almost singular" it is "almost noninvertible." This means 
that there are difficulties in finding the contravariant coefficients from the covariant coefficients, since this requires 
the inverse of G n k- The implication of this is possibly of pragmatic importance: the near singularity makes it difficult 
to find the "theoretical" coefficients from the "experimental" ones. 



III. INDUCED COMPLETENESS 



Quantifying QN excitations 



8 



A. Small change 

Our approach to quantification of excitation is based on the idea of complete QN systems. It is important to 
demonstrate that complete QN systems exist and to have an example of a system in which we can compute excitation 
coefficients. We must also deal with a more specific issue. The QN modes of potentials for black hole problems do not 
form a complete set; both the Regge- Wheeler potential, which describes perturbations of a Schwarzschild black hole, 
and the Poschl- Teller potential [pLlf , which has been used as an approximation to the black hole potential have 
quasinormal mode sets which cannot completely describe the wave form resulting from an initial perturbation. An 
important question is what the relevance is to black hole processes, of any quantification based on an assumption of 
complete systems. The details of this section are somewhat disjoint from the main goal of the paper: quantification 
of QN excitation, but are a necessary step in developing our approach. (They are also rather interesting in their own 
right!) 

The key idea in developing a demonstrably complete QN system, and in showing its relevance to black holes, is 
"induced completeness." It has been argued [Q that for the problem defined by (||), the eigenmodes will be complete 
if there are two values of x at which V(x) is not C°°. If this is true, it appears that the black hole potentials can 
be modified in such a way that completeness of QN oscillations is induced, while the effect on other aspects of the 
problem - in particular on the evolution of Cauchy data - is kept arbitrarily small. In this section we explore induced 
completeness. 

It is not our purpose here to look into the generality of induced completeness. Rather, we want to have a specific 
example of a complete QN system with which to explore the question of quantifiability of QN excitations. We will 
therefore focus on a very specific model that contains most of the flavor of black hole potentials, but is simple enough 
to allow rather straightforward analysis and a proof of completeness. 

We use the fact that (||) can be solved with elementary functions in the case that the potential has the form of a 
truncated centrifugal potential, 

V TCp(x) = { e{£+1)/x 2 I til ' 

where I is an integer. It gives a good representation of the sharp cutoff of the black hole potential as x — > — oo and 
the approximate asymptotic form of the black hole potential at x — > oo. It differs from the true black hole potential 
in the details of the x — > oo potential that give rise to the power-law late time tails of black hole waveforms. These 
tails almost certainly are an obstacle to an analysis of QN excitation, and a modification of the potential to eliminate 
these tails has already been made in work on QN excitation. 

We will hereafter consider only the case I = 1, to be called the "truncated dipole potential" (TDP). An example 
of the simplicity of the TDP problem is that it has only a single pair of QN frequencies: ui\ = (1 + i)/2xo and 
ujf = (— 1 + i)/2xo (see Appendix IA). There are several ways in which we could try to induce completeness into this 
problem. We have studied both the addition of a small step (with discontinuities) to the TDP and the addition of 
a delta function. The delta function has the disadvantage of its distributional nature, but it offers the advantage of 
considerable simplicity as compared with the step. We have found no significant "practical" difference between the 
results (QN locations, convergence of QN sums) between the two examples, so we will describe here completeness 
induced with a delta function. The total potential in this case will be called the "spiked truncated dipole potential" 
(STDP) and is given by 

VvrD P (x) = VTDr(x)(l + V B 8(x-xs)). (18) 

where Vtdp (x) is the £ = 1 form of the potential in ([lT]) . 

We first establish that the influence on the evolution of initial data can be made arbitrarily small. To do this we 
choose the Cauchy data, at t = to be given by 

V>o = sin ( 2?r - - X2 ) (x 2 < x' < x x ) (19) 
V Xi-x 2 J 

ipo = -dip /dx , 

with xi < xo, so that the initial data has the form of one full cycle of a sine wave, located to the left of the potential, 
traveling to the right. (This Cauchy data is chosen for convenience in demonstrating mathematical points; it has no 
justification as initial data for gravitational waves being produced in the neighborhood of a black hole.) 

For the results to be shown here and in following sections, we choose Xq = 1, x% = 1, x 2 = —5, and x$ = 10, and 
^obs = 120. For this choice of X\ and x 2 , the sine has a wavelength, 6, that is roughly half that of the QN oscillation 
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of the TDP . This allows us to see similar, but distinguishable signs in the evolved waveform of the propagation of 
the Cauchy data and of QN oscillation. 




FIG. 4. Time evolution for initial data incident on the TDP and on the TDP with an added <5-mnction. Left: linear plot, 
right: logarithmic plot 

Figure || shows the time evolution of the Cauchy data for the original TDP and for the TDP with an added 5- 
function with different amplitudes ranging from Vg = 1 to Vg = 10~ 6 . The waveforms are followed out to times at 
which they have decreased in magnitude from the maximum by a factor of ~ I0~ 15 , at which point numerical error 
obscures the results. Even for the largest value of Vg, the influence of the delta function on evolution is visible in the 
linear plot only after we change the scale for f(t) at t — 131. To see the effects more clearly we also plot the logarithm 
of \f(t)\ for a longer observation time, showing several interesting features. For all but the largest amplitude of the 
delta function (Vg = 1) the waveform after the first few oscillations consists only of QN ringing, with the characteristic 
distance At = |3fiw7|7r = 2ttx — 2tt between zeroes of f(t). For Vg — 10~ 3 , the effect of the delta function shows up 
only as a phase shift after the the first four, or so, full cycles of oscillation. The effect of the smallest delta function 
amplitude Vg — 10~ 6 is smaller by two orders of magnitude than that for Vg = 10 -3 , and too small to be seen even 
in the logarithmic plot of Fig.|]. 

These results make it clear that any reasonable measure of the influence of the delta function, such as the integral 
of the square deviation from the TDP waveform, is tiny and decreases with decreasing Vg. We will accept these 
numerical results as a sufficient demonstration of this point, and will not attempt an analytic proof of this point. 



B. QN spectra 

Although the influence of a small delta function on the evolution of Cauchy data is negligible, the influence on the 
spectrum of QN frequencies is profound. The method of computing the QN frequencies for the STDP is outlined in 
Appendix I. The results of the computation are presented in Fig. || for the same potentials considered in Fig. |J 

The original TDP has only one pair of QN frequencies at (±1 + i)/2. With an added 5-function, an entirely new 
set of modes appear. Note that the real parts of the additional frequencies seem to be unbounded. The asymptotic 
spacing 3?(u; n ) = 27r/10 of the real parts of the frequencies is shown in Appendix IC to correspond to the length of 
the "cavity" bounded by Xq and the 5-function: L = xg — x — 10. The imaginary parts are shown in the Appendix 
to increase as the logarithm of the real part. 

When Vg becomes sufficiently small, one of the QN frequencies approaches the value (l + i)/2xo of the "native" QN 
frequency, that of the original, TDP . This can be seen more clearly in Fig. which shows the "path" of this QN 
frequency in the complex plane for values of Vg varying in steps of 1/2 from Vg = 1 to Vg ~ 10 -6 . As Vg decreases, 
the imaginary parts of the additional frequencies increase, moving them away from the native one, eventually leading 
to two very distinct subsets of QN frequencies. 
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FIG. 5. Quasinormal frequencies for the TDP with an ad- FIG. 6. "Path" of the original quasinormal frequency of 



ditional 5-function at x& = 10. 



the TDP as the amplitude of the additional 5-function is de- 
creased. 



It might seem that it is an obvious necessity for the QN spectrum to have a mode approximately at the location of 
the native mode, since the evolution of Cauchy data is only slightly affected. This turns out not to be true, however, 
for other ways of inducing apparent completeness |]i~3f . Cutting off the TDP potential at some very large value of x, for 
example, has a negligible effect on the evolution of Cauchy data, and it also produces an unbounded set of additional 
QN frequencies. However, the spectrum of QN frequencies turns out to contain no frequency near the location of the 
native mode. 



C. Complete sums 

We now develop the connection between the Cauchy data tpo and ipo and the coefficients of QN oscillations (i.e., the 
contravariant components a n in the case of complete QN bases). Here we simply outline how coefficients associated 
with QN basis functions are computed in general. In Appendix II a proof of completeness of sums with these coefficients 
will be given for the case of the spiked TDP, and of Cauchy data meeting certain criteria. One of the criteria will be 
compact support for the Cauchy data, and we will assume from the outset that i/'o and ipo vanish outside the interval 
Xi < X < x±. 

The QN coefficients are found by starting with a Laplace transform 

POO 

ip(s,x) = / e~ st V(t,x)dt . (20) 
Jo 

(The connection between the Laplace variable s and the Fourier ui used in Appendix I and most of the paper, is 
through the correspondence s ^ iui.) The transformed wave equation reads 

di>{s, x)/dx 2 - [s 2 + V(x)] 4>(s, x) = J(s, x) , (21) 

in which the source J(s, x) is determined by the Cauchy data: 

J(a,x) = -i>o{x) -sVo(z) (22) 

A solution can be found in the form 



ip(s,x) = / G(s,x,x')J(s,x')dx' 

J X 2 



(23) 
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Here the Green function can be constructed from the homogeneous solutions /_(s, x) and /_(s,x) of of the wave 
equation (f2l|): 



G(s, £, x') 



f-(s,x')f + (s,x) (x' <x), 
W(s) {f_(s,x)f + (s,x>) (x'>x), 



(24) 



where W(s), the Wronskian of /_ and /+ is independent of x and x' . 

We will assume that the potential falls off quickly enough at x — » ±00 so that homogeneous solutions /_ and /+ 
can be found with the property 



f-(s,x) ~e sx (l + 0(l/\x\)) as x 



and 



f+{s,x) 



(l + 0(l/|x|)) as x^oo 



(25) 



This condition is satisfied for black hole potentials and for the TDP. 

Once we have found a solution of the Laplace transformed wave equation, we can reconstruct a solution of the 
time-dependent wave equation by applying the inverse Laplace transform: 



^(t,x) = — ( e st ?p(s,x)ds 

where F denotes the path of integration, which lies parallel to and just to the right of the imaginary axis. 

Im(s) 



(26) 




FIG. 7. Closing the contour for integration in the complex s plane. The curve F is the original contour of integration. It 
can be closed by the addition of an arc at infinity. If there is a branch point at the origin, as shown, then a branch cut can be 
drawn from s — to s — —00, and the contour shown can be drawn. 

Figure^ illustrates features in the complex s plane. There may be poles and essential singularities, drawn as xs, 
and there may be a branch point, like that shown at s = in the figure. The contour T can be closed with the addition 
of a single arc to the left if ip(s,x) is a single valued function of s to the left of T. If there are branch points then 
a path like that shown in Fig.^ must be used. The Green's function for the Regge- Wheeler potential has a branch 
point at s = 0, while the TDP and the STDP do not. 

The integral for ip may then be evaluated as 



$(t,x) = V'Res(e s *i/)(s,a;),s : ,-) + — / e st 4>(s, x)ds 

2tti J c 



(27) 



where the first term is the sum of residues at the singularities inside the completed contour. The second term represents 
contributions to the integral from arcs at 00, and along branch cuts. In the case of the spiked TDP there is no branch 
cut, and the arc at 00 makes no contributions, so we are left only with the first term, a sum of oscillations at discrete 
frequencies Sj that correspond to singularities of %jj{s,x). From ( |23| ) and (|24|) we see that singularities in i>{s, x) must 
either be singularities of the homogeneous solutions /+ , /_ or zeroes of the Wronskian. For any finite s wc can find 
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homogeneous solutions /+ and /_ of (21); therefore, singularities in the Green's function can occur only at zeros of 
W(s), and the residues in (5?]) must be taken at these zeros. But the vanishing of W(s) at Sj means that / + oc /_ at 
So , and hence that s, is a QN frequency g. 

If, as in the case of the sTDP, the only contributions to (^?]) occur in the sum, then we are left with 

■qi(t, x) = ^ Res(e s *^(s, x),Sj) , (28) 

where the sum is over QN frequencies. 

We now assume that the zeros of W(s) are all first order, so that 



urr . dW 
W{s) = — 

as 



( s - S] ) + 0[(s-s 3 ) 2 }. (29) 

Sj 

For x > x\ (i.e., for x to the right of the region in which the Cauchy data have support) we have, from ([24] ) that 

G(s,x,x') = L^lM+i^ +ois-sj) (30) 



and hence, from (E3j), the residues of -0 at s ? is 



Res(^(s,x) = ^J^k¥^ I J{sj, x')f-( Sj , x') dx' . (31) 



dW/ds\ S] 

Finally, the sum over the QN basis functions can be written 

*(t,z) = £> % (*,z) , (32) 

where 

Uj (t,x) = f+{s ] ,x)e^ t (33) 



and 

1 



dW/ds\ 



(jj {x') + SjtPo^')) U(sj,x')dx' . (34) 



The waveform at x b s > x\ is then given by (|) with 

iuj-i = Sj a 3 = b> f + (sj,x obs ) . (35) 

This prescription has been applied to the spiked TDP, and the sine wave Cauchy data, with our standard choice of 
parameter values, Xq = 1, xg = 10. Notice that our choice a; 2 = —5 in jl9| ) satisfies the criterion in the convergence 
proof that the continuous, but nonsmooth, Cauchy data have support only for x > xq — 2/3(x$ — xq) = —5. For this 
potential and these Cauchy data, the functions rpo(x),ipo(x) and f-{sj, x) are trigonometric or exponential functions, 
so the integral in (|34|) is elementary. Once the QN frequencies, u>j, and the factors dW/ds\ alj have been computed 
(as described in Appendix I), the coefficients are easily evaluated. The figures below show the result of using these 
coefficients in sums of the form (||). 

Figure [| shows the computed result for time evolution of the Cauchy data in the case that V§ = 1. This is compared 
along with mode sums for an increasing number of terms, up to N = 10, 001 terms. Figure^] shows the same plots in 
the case that Vg = 10 -6 . In order to avoid cluttering the picture, we do not plot the values of the mode sums when 
they are far from having converged. These figures illustrate the fact that the mode sum converge faster for a larger 
(5-function. Also, for the smaller (5-function, the mode sum converge more rapidly at later times, with convergence 
"sweeping down" from late to early times. 
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Time evolution 
and N > 101 




120 125 130 t 120 125 130 t 

FIG. 8. Values of the mode sum for different number of FIG. 9. Values of the mode sum for different number of 
terms N, compared to the wave form resulting from integrat- terms N, compared to the wave form resulting from integrat- 
ing the time dependent wave equation (Vs = 1). ing the time dependent wave equation (Vs = 10 -6 ). 

For a better view of the differences between the mode sums and the time evolved wave functions, Fig. ^ shows the 
logarithmic difference (evolved waveform vs. mode sum) for Vs = 10~ 6 . 

Some of the systematics shown in these figures can be heuristically understood. The convergence in the Vs = 1 case 
is similar to that of a fourier series (equally rapid at different times) since the QN frequencies, in this case, with small 
imaginary parts, are similar to the real frequencies of a fourier series. On the other hand, for Vs = 10 -6 , the damping 
of the additional (i.e., non-native) modes is much stronger and increases faster with N. Any error that originates 
from cutting off the mode sum after a finite number of terms can be regarded as being composed of modes with very 
strong damping; these modes are very large at early time. 

We point out here an interesting technical feature of these results. The mode sums are so accurate at later times that 
the differences shown in Fig.[l(] for N = 10,001, is actually dominated, after t ~ 122, by the numerical truncation 
error in computing the evolved waveform. A smaller time step in numerical evolution can improve the numerical 
accuracy of the computed waveform, and can move to a slightly larger time the point at which the evolution vs. sum 
difference is dominated by the numerical errors in the evolution. 



/(*) 




120 125 130 t 



FIG. 10. Logarithmic difference A/ between the values of 
the mode sum for different numbers of terms and the wave form 
resulting from integrating the time dependent wave equation 
(Vi = lO" 6 ). 
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IV. QUANTIFYING EXCITATION 
A. The excitation coefficient 

We now return to the question of how to quantify the excitation of QN oscillation. In Sec. II we defined the 
excitation coefficient Ak for a QN oscillation in a complete QN system. Due to the time shift problem, we argued that 
the excitation coefficient seems the only plausible indicator of the QN content of a waveform. In Sec. Ill we have seen 
that, at least in the particular example of the TDP, completeness can be induced in an incomplete system to create 
one that is "physically equivalent," i.e., differs negligibly in the evolution of Cauchy data. We can now ask whether, 
at least for the model problem at hand, we can use the excitation coefficients of the completed system to quantify the 
QN oscillation in the original ("native") system. 

It is worth emphasizing that the TDP is a particularly simple starting point for these considerations, since it 
has only a single conjugate pair of native QN modes. The sine wave data we have used in the previous sections is 
also convenient since it produces a waveform (see Fig.^J) which clearly contains QN ringing, but contains significant 
oscillation at a different frequency (u = 2w/6). 

The results of computations with this model are presented in Tabic |. Results are shown for the spiked TDP for 
different values of Vg and for the standard Cauchy data, a right moving sine wave initially extending over the interval 
[—5,1]. For comparison, a shifted sine wave, initially at [—8,-2], was also computed in the case Vs — 10 -6 . Also 
included are results for the smooth TDP, and for the Zerilli equation with the Regge- Wheeler potential with initial 
data corresponding the the Close Limit technique for black hole collisions jl4j . 

The table presents values of excitation coefficients A\ = a 1 (ai)* + (a 1 )*(ai) and the energy excitation coefficients 
Ei = (wj) 2 (ai)*a 1 + (wi) 2 (aj)*a 1 , introduced in ( [To]) and ([15]). In each case the QN frequency uj\ is taken to be 
that eigenfrequency which corresponds to the native mode in the limit Vs — > 0; that is, io\ always lies on the path 
in the complex plane shown in Fig.[(| For small values of Vs this QN frequency is close to the native QN frequency 
(±1 + i)/2. The waveform norm and the total energy are computed from wave forms obtained by explicit numerical 
integration of the wave equation (|l|). Also included are estimates for the numerical uncertainties of the results. 

The norm and total energy can be computed from Eq.(|^), or directly from the wave form. Similarly, the covariant 
coefficient may be computed from the contravariant ones, using Eq.(fij), or from the wave form via Eq.(@). Even the 
contravariant coefficients can be obtained from the wave form itself, using an asymptotic fit as it approaches a pure 
QN oscillation at the least damped QN frequency, instead of the residues of the Green's function through Eqs. ( |34| ) 
and (|35|). The wave form itself, in turn, can be obtained by direct numerical integration of the wave equation, or by 
using the mode sum as in Eq.(||). We have checked these alternatives for obtaining the quantities listed in the table; 
they result in essentially the same values as the route we have taken here. 

We first notice that the values for the excitation coefficients are identical, within the numerical errors, for the 
original initial data (sine wave over the interval [—5, 1]) and the shifted initial data (over the interval [—8, —2]). This 
confirms our earlier argument that the excitation coefficients we have defined are independent of a translation of the 
initial data or, correspondingly, to a time shift of the evolved wave form. 



V(x) 


V s 


Initial data 


Ai 


E x 




/ 


2 




/ 


2 


STDP 


10" a 


Sine [-5, 1] 


-40 854. ± 1. 


-23 810. ± 1. 


19 246. ± 7. 


27 404. ± 12. 


STDP 


1CT 6 


Sine [-5, 1] 


-57 933. ±1. 


5 804.7 ±0.1 


19 246. ± 7. 


27 404. ± 12. 


STDP 


1CT 6 


Sine [-8, -2] 


-57 934. ± 1. 


5 804.5 ± 0.3 


19 246. ± 7. 


27 404. ± 12. 


STDP 


io- 7 


Sine [-5, 1] 


-57 879. ± 1. 


5 811.9 ±0.1 














STDP 


10~ 8 


Sine [-5, 1] 


-57 873. ± 1. 


5 812.6 ± 0.1 














STDP 


10" 9 


Sine [-5, 1] 


-57 873. ± 1. 


5 812.6 ± 0.1 














TDP 




Sine [-5, 1] 


-57 872.8... 


5 812.64. . . 


19 248.0 . . . 


27 415.5... 


Zerilli 




Close Limit 


10.699 x 10~ 4 




5.566 x 10 -4 





TABLE I. Excitation coefficients and norms for the spiked TDP QN system and for black holes in the close limit approxi- 
mation. 
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It is important to realize that there are different ways to compute the quantities in the table, and several of them 
do not depend on the complete set of modes. The waveform and energy norm, of course, require only the waveform, 
but all coefficients can also be computed entirely from the waveform itself. The least damped QN frequency can be 
inferred from the asymptotic late time behavior, as can the contravariant coefficient for the least damped mode. [See 
(0).] Of the cases we study the least damped mode always corresponds to the native mode except for the Vg = 1 
model. (See Fig.||.) Once the QN frequency is known the covariant coefficient can be computed directly from the 
integral in (Q) and, with the contravariant coefficient known, the excitation and energy coefficients can then be found 
from ( |l3| ) and (^5|). Since the quantities in the table can be computed from waveforms, we can compute them for 
the smooth TDP. In this case we do not have a complete set of modes, but that is irrelevant to the procedure for 
computation. (It turns out, in fact, that the simplicity of the smooth TDP and the sine wave data allows a closed 
form solution for the waveform, and for the norms and coefficients. This closed form solution has been used, and the 
values for the smooth TDP can be found for the table to essentially arbitrary precision.) 

The numbers in the table make it clear that the results for the smooth TDP are the Vg — ► limit of the spiked 
TDP. This is obvious from a computational point of view, since all results can be computed from the waveform, 
and the Vg — ► waveform approaches the smooth TDP waveform. From another point of view, however, this result 
is interesting and important. It means that we can compute the excitation coefficient independent of the method 
with which completeness is induced. Put another way, it means that we can compute the excitation coefficient for a 
small completeness-inducing change, independent of the nature of the change. This conclusion, in fact, is crucial to 
the possible importance of the excitation coefficient. We could not use the excitation coefficient to characterize an 
excitation of a physical system if the value of the coefficient depended on our choice of a modification of that physical 
system. 

The table also gives a value for the excitation coefficient for the gravitational radiation produced, in the close limit 
approximation [ fijj , by the head on collision of two equal mass nonrotating holes. In this case the QN spectrum is 
infinite, but the QN oscillations are not complete. The excitation coefficient given is for the least damped of the QN 
oscillations, a frequency that appears to dominate the waveform produced. 

Up to this point we have noted that computations with the model problem illustrate and confirm the features of the 
excitation coefficient that made it an attractive candidate for the quantification of the excitation of QN oscillations. 
Table|, however, also makes it unlikely that the excitation is a useful index of QN oscillation. For the small-(5 spiked 
TDP the excitation coefficient is negative, and is larger (by a factor ~ 3) in magnitude than the norm. Note that 
we cannot ignore this "wrong" sign, and simply keep the large magnitude as an indication of a strong QN presence. 
Due to relation (^) we must conclude that the sum of all the other QN oscillations (those unrelated to the native QN 
mode) must be greater (by a factor ~ 4) than the norm. For the close limit waveform the excitation coefficient is 
roughly twice the size of the norm of / and hence the sum of excitation coefficients of all other QN oscillations (if the 
system were made complete somehow) would be negative. 

As a possible alternative to the excitation coefficient we have also computed the energy excitation coefficient, as 
defined in (jlj). The results listed for E\ and the norm of /, however, do not make this any more attractive as a 
measure of QN excitation. 

At this point, we also note that the excitation coefficient A\ is not related to the norm of the QN mode contribution 
corresponding to the first pair of modes, i.e. of b x u\ + {b x u\)* . Similarly, the energy excitation coefficient E\ is no 
measure for the energy of this QN mode contribution. 

It is, of course, impossible to prove that something like quantifying QN excitations cannot be done in a mathemat- 
ically acceptable way. Despite this, and perhaps to provoke further work (by others!) we are tempted to offer the 
following very tentative conclusion: Consider the following three conditions, which allow a mathematically meaningful 
as well as a useful measure of excitation; they are all satisfied for normal mode systems: 

• The measure of excitation is independent of a simple shift of the wave form (i.e. a shift in space of the initial 
data, corresponding to a shift in time for the time evolved wave form) . 

• Excitation strength can be quantified individually for any number of modes, with the individual excitations 
adding up to the total norm of the wave form. 

• The measure of excitation is useful for quantifying the excitation in a comparative way; in particular, it always 
lies between 0% and 100%. 

We conjecture: There is no quantification of QN oscillations of the waveform, based on the algebra of the function 
space of QN oscillations, which satisfies all three of these criteria. 

The excitation coefficients we are defining in this paper satisfy the first two conditions, but not the third one (In 
Table| we see that the excitation coefficient can be negative and can be greater than 100%.) One might regard these 
two conditions as related to mathematical properties of the QN mode system, while the third one is of a more practical 
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significance. We are currently investigating another technique to quantify excitations, with a measure that satisfies 
the first and third criteria, but not the second one. This measure may turn out to be of some practical utility, but 
is not as closely related to the mathematics of the function space as are the present considerations. A description of 
this work will be published elsewhere. 

Our conjecture is based specifically on the appearance of the excitation coefficient as the only quantity in that 
function space that solves the "time shift" problem, and on the observed failure of the excitation coefficient to be 
"useful." It is also based, less specifically, on our belief that there are differences between normal mode and QN 
systems that cannot be bridged. For this reason we speculate that inducing completeness, while it is mathematically 
interesting, is not likely to lead to useful tools for understanding the underlying native system. The reason is that 
the spectrum of added QN frequencies is unrelated to the native system and characterizes only the method used to 
induce the completeness. We, in fact, are willing to extrapolate in the direction of this speculation. We conjecture 
that in a complete QN spectrum which has a dense set of frequencies, and a small number of "isolated" frequencies, 
like the spectrum of the spiked TDP s in Fig.|B|, it is useful to modify the system to remove completeness, in order to 
get a more useful understanding and a simplified method of analysis. The obvious example of this is removing the S 
functions from the spiked TDP problem in order to get a physically equivalent, but simpler incomplete system. 

A very similar point of view is that in a QN spectrum, not all frequencies are equally important. Some will actually 
be evident in waveforms produced in the evolution of generic Cauchy data; others won't. For the small-<5 spectra of 
Fig. [|, the "interesting" QN frequencies are the isolated ones near (1 +i)/2. In the Vg = 1 case, on the other hand, there 
is again the appearance of QN ringing in the waveform, but the spectrum contains no isolated frequency. A similar 
situation was found in a study where the Regge- Wheeler potential was replaced by a series of step potentials 0|. 
More generally, one should ask: If one has only the spectrum and the associated quantities (e.g., the metric matrix), 
is there a way of identifying which QN frequencies are "important" in the sense of really characterizing the evolution 
of Cauchy data? In this sense we are asking a question related to "to what extent are (some) QN modes like normal 
modes?" since normal modes do characterize the system in which they arise. 



B. Condition number of the metric matrix 



The metric matrix (|TTj) would appear to be a likely place to find a way of characterizing a QN system without regard 
to specific waveforms produced by specific Cauchy data. We will now discuss some numerical res ults which confirm 



the intuitive insight we had gained by doing thought experiments on specific Cauchy data in Sec. IIP. Studying the 



metric matrix directly, we don't need to refer to specific initial data any more, as we had to do before. 

We first note the singular nature of the infinite metric matrix in two cases. To characterize this singular property we 
truncate the set of QN functions, keeping only the first D. We then compute the condition number R (ratio of maxi- 
mum to minimum eigenvalue) for the D dimensional subspace. The condition number as a function of -D S ub is plotted in 




2 4 8 10 20 D suh 



FIG. 11. The condition number R for the metric matrix for 
a truncated subspace of dimension D su b, spanned by the first 
D su b QN basis functions. Results are shown for both the spiked 
TDP (Vs = 1(T 6 ) and the Regge- Wheeler potential. 



Fig. |TT] for two cases: the spiked TDP and the QN spectrum of Schwarzschild black holes (the QN modes of the Zerilli 
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or Regge- Wheeler potentials). For the spiked TDP the approximately straight line in the log-log plot suggests that 
the condition number increases roughly as the twelfth power of the dimension -D su b of the subspace. For the black hole 
QN spectrum the increase is even more dramatic, suggesting perhaps a "more singular" metric for this incomplete 
QN spectrum. 



C. Angles between basis functions 



In Sec. II we introduced covariant basis functions <fi m (t) with the definition 4> m (t)-e 



iu) n t 



If the basis functions 



were orthogonal we would have that <fr n (t) and e luJnt are "parallel," that is, <j) n (t) oc e 1 "™*. It is plausible that 
such statements as "this wave is dominated by oscillation at the fundamental QN frequency" are most meaningful if 
the covariant and contravariant basis vectors are nearly "parallel." To measure the extent to which 4>>(t) and e luJjt 
fail to be proportional we can introduce an angle ctj between them, defined by 



cos(aj) 



1 



(36) 



where G lJ is the matrix inverse of G; 
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FIG. 12. Cosine of angles between several contravariant FIG. 13. Cosine of angles between several contravariant 
and covariant basis vectors, as a function of the dimension of and covariant basis vectors, as a function of the dimension of 
a subspace spanned by D su b QN modes of the spiked TDP a subspace spanned by Asub QN modes of the Regge- Wheeler 
with Vs — 10 -6 . potential. 



The components of the infinite matrix G lJ cannot be computed, so again we truncate a subspace by keeping only the 
first D S ub vectors, and we compute the angles in that subspace with (^6|). Results are shown in Fig. H| for the spiked 
TDP with Vs = 10~ 6 . The value of cos (a, ) is shown for several QN modes as a function of -P S ub, the dimension at 
which the subspace is truncated. Figure [l3[ shows cos (a,) for the black hole QN spectrum. In Fig. [IJ, the decrease of 
cos(ai), corresponding to the native mode, is much slower than that of the additional modes. One might speculate 
that this is related to the fact that the native mode is more characteristic of the system than the additional ones. 
No such clear distinction can be seen in Fig. [p| for the case of the Regge- Wheeler potential. However, the angles 
increase much faster for the more highly damped modes as well, which might again indicate that the fundamental, 
least damped mode is more characteristic of the system than the more strongly damped ones. 
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APPENDIX I. FINDING QUASINORMAL FREQUENCIES OF THE TDP AND OF THE SPIKED TDP 

A. The unmodified TDP 

The TDP is defined as 

V W = \l(l + l)/x 2 x>x ' ( 3? ) 

where Z is an integer. We look only at I = 1, but the procedure can easily be extended for larger values of I. Also, we 
will generally let xq = 1. 

The domain of the wave equation is naturally divided into two regions: 

I. x < Xq In this region, the potential vanishes, and therefore the wave equation has the two trivial solutions: 

^n(x) = e +iux (38) 
M x ) = z~ lUX (39) 

II. x > xq For integer values of I, the solutions are given by finite sums. For I = 1, we have 

Vm(x) = e+^(l - — ) (40) 
iujx 

ipm(x) (1 + — ) . (41) 

IUJX 

Obviously, the solutions satisfying the required boundary conditions at negative and positive infinity are 

ip-(x) = tpn(x) {x < x ) (42) 
ip+(x) = ipm(x)(x >x Q ). (43) 

In general, of course, ip-( x ) wu l be a linear combination of ipm(x) and i/^uaix) for x > xq, and ip+( x ) a combination 
of ipn(x) and tpi2(x) for x < xq. A quasinormal mode is a solution where both boundary conditions are satisfied 
simultaneously, i.e. ip-( x ) — i>+{ x )- The easiest way to find out if this is the case is to compare ip-(x) and ip-\-(x), 
as defined in (pl2|), at x — xq. Strictly speaking, ip-( x ) is n °t defined at x = xq. However, any solution of the wave 
equation will have to be continuous, and have a continuous first derivative, at x = Xo. It is therefore permissible 
to use the left limit of ')p-(x) and of ip'_{x) as x — > Xq- , and compare them with the values of ip+( x ) an d a ^ 

x = Xq. 

The comparison is done using the Wronskian determinant of rft_(x) and ij)+(x): 

W[ip-,ip+}{x ) = i(;-(xo)ip' + (xo) - ip'_(x Q )^ + (x ) 

2 1 

= -2iu «. (44) 

xo iujXq 

Solving the equation W[tp- 1 tp+](xo) — for u> yields the quasinormal frequencies 

±l + i 

D = 

2x 

Therefore, for 1 = 1, there is only one pair of quasinormal frequencies. 



(45) 



B. QN frequencies of the spiked TDP 



We define the 'spiked' potential as 

V(x) = V(x) + V s S(x - x s ) 

We now have to distinguish three areas: 



(46) 
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I. x < xq This region is identical to region I for the unmodified potential, with the same set of two solutions. 
Ila. xq < x < Xg Again, there are two linearly independent solutions 

^iiai(z) = e +i " x (l - — ) (47) 

IIVX 

ifartix) (1 + — ) . (48) 

ILUX 

lib. x > x§ The two linearly independent solutions are 

ifobi(x)=e* u *(l-^-) (49) 
uvx 

ipub2(x) = e- iux (1 + — ) . (50) 

IU)X 

In our notation the function ijjuh2(x), for example, refers to the solution in all regions that, in region lib, has the 
functional form shown in Eq. |5(]. 

Due to the S function separating regions Ila and lib, the solutions ipiia,2( x ) and ipub2(%) are not the same. Rather, 
■0iib2(^) will be a linear combination of ipua.i( x ) and ipu^ix) in region Ila: 

1pllh2(x) =PlVTlal(a;) + P2lpIIa,2(x) (51) 

Once again, the solutions 

%l)-{x) = ipn{x)(x <x ) (52) 

ip+(x) = ipiib2(x) (x>x 5 ) (53) 

satisfy the boundary conditions at negative and positive infinity. However, in order to compare solutions at x — xq, 
we now have to determine the representation of ipiib2(x) in region Ila, i.e. we need to know the coefficients for the 
linear combination (fill). 

These coefficients can be determined using the junction conditions for any solution ip(x) of the wave equation across 
the S function: 

1. ip(x) must be continuous at x$ 

2. the derivative ip'(x) must have a discontinuity given by 

iP'(x s +) - = Wfo) (54) 

The second condition is obtained by integrating the wave equation from x = x$ — e to x = x$ + e, letting e — > and 
using the first condition. 

Using ipub2 for ip'(x s + ) and Eq. ( |5l"|) for ip'(xg~), we can solve these conditions for pi and p 2 . We obtain 

T/- [^Ha2(%)] 2 

Pi = Vs — TTr (55) 

vv+12 

p2 = 1 _yMf^™(*s) (56) 
W + i2 

where W +12 = W[^i B i, ihiaa] = 
Therefore, 

W[tp-, ip+](x ) = W / [V'n,i/'iib2](a;o) = PiW[ip n , ^HaiR^o) + P2W[ipii,ihua\(x ) 

= Ri + V s (R 2 + R 3 e- 2tuL ) , (57) 

where L — x s — xq, and W[tpn, V'liaiK^o) = e 2lulXo / '(iujx 2 ,), and thus 
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Ri = W[ipn,ipu a2 ]{xo) = -2ilu - — - (58) 

i? 2 = - ^axC^^K) ^ ^ = - f 1 - 7^) fl + A- + (59) 
VK+12 \ [iwx 5 y J \ iujxq 2{iux yj 

Rs = e 2 - i %^W I1 ^n al ]K) - -^T" 1 ^ f 1 + t^-Y (60) 

2(iu}x y \ iujx s J 

The quasinormal frequencies of the spiked TDPcan now be computed numerically by searching for roots of the 
equation 

W^_,V+](x o )=0 . (61) 



C. Asymptotic approximation for the QN frequencies of the spiked TDP 

It is possible to find an asymptotic formula for the QN frequencies under the assumption that the absolute value 
of the frequency becomes large. We start by assuming that in (|6l]), we have \uixo\ ^> 1 and \u>xg\ 3> 1. The condition 
for QN frequencies can then be written 



2iLux + 0([ux)} ) + V s 



1 + OHlux)]- 1 ) + ^7—32 (1 + OHujx)]- 1 )) 



= , (62) 



2(iwxo) 

where x is the minimum of xq and x$. It is clear that for \u>x\ ^> 1 there can be solutions only if 

and we use this approximation to find an iterative solution for the QN frequencies. We start by taking the cube root 
of 

e -2i«L = 4 M^0 (64) 

V S 1 ' 



to write 



e -|fa,«£ e imL e *A = _ f^o\ 3 lL0 = _A iuj , (65) 



Vs 



where A = j§7r, j = 0, 1, 2. 

Taking the absolute values on both sides gives 



e^' L = A\u\ (66) 



With u> — lur + iu>i this last relation already shows that uii <C ~ lor is required for a QN frequency. 
Using ( |66| ) to rewrite ( |65|) we find 

cos(§Lw fl + A) - i sin(§Lw fl + A) = -i^- (67) 

O O \LO\ 

We now make the approximation ujj ks 0, i.e. oj « cj^ |cj|, leading to 

2 2 
cos(-Lwr + A) = 0, sin(-£wR + A) = 1 (68) 

= | + 2J\T7r - |j7r, N = 0,1,...; j -0,1,2 (69) 

or equivalently 
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1 37T 

{io R ) n = -{— +mt), 71 = 0,1,... (70) 

An approximation for ujj is then obtained by using (|(36|): 

3 1 

The approximate solutions of ( |70|) and (|7l]) can now be iteratively improved. We rewrite ([n]) as 



V S R 3 



and take absolute values to get 



2uj r L 



R(w)=R(w R ,uj) (72) 



\R\ (73) 



With these inserted in IY72) we arrive at 



" I= 2L (74) 



COS (2w R L) = -^0-, sin(2wfl L) = ^. (75) 
This can be used as an iterative method to find the p th iteration from the p — 1 th approximation, as follows: 

cos(2(ujr) f L) = ; — — — z ~ — r — t— — , sin(2(ujjiyL) = — — ; — — (76) 



and 



= ln|fl((o; fl )P,(a; J )P- 1 )| 
1 I] 2L 



(77) 



This iterative solution can be started with any value of n in the zeroth approximation of (|7^) and (|FH), and the 
iteration converges to the exact solution. The iteration cannot be used to find the native QN frequency of the smooth 
TDP since V$ = in that case. 

APPENDIX II. PROOF OF CONVERGENCE 

A proof is given here of the convergence of the sum of quasinormal excitations for the spiked TDP under approprite 
restrictions on the Cauchy data. To do this we start by defining the integral 

f 1 f s 2 d>( s)e s( - dl+d ^ 

*w a / *.) ds - (78) 

Here A(s) is defined to be 

A(s) =P 1 (s)+e 2sL s 3 P 2 (s) , (79) 

in which Pi and P2 are polynomials of finite order in 1/s: 

P 1 {s) = A x +B 1 /s + --- P 2 (s) = A 2 +B 2 /s + --- . (80) 

The path of integration in the complex s plane is along the vertical axis, from — ioo to +ioo. For 3?(s) < the function 
<fi(s) is required to have the property that for some real constants and p 
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(This condition will be related below to restrictions on acceptable Cauchy data.) The constants appearing in J78| ) 
are taken to satisfy the following conditions: (i) The ratio Ai/A 2 is real and positive, (ii) L,d\ and d 2 are real and 
nonnegative, and 2L > d%. (iii) The constant p appearing in must be large enough that 



M2 
2L 



2 > 



(82) 



The roots of A are denoted Sk- (They represent, of course, the QN frequencies according to the usual correspondence 
s <-> iuj.) . Since they must occur in complex conjugate pairs it is convenient here to use the notation s±, s_i, S2, s_2, ■ ■ ■ 
with s±k indicating the corresponding roots with positive and negative imaginary parts. 




Im s 



Re S 



FIG. 14. Contour for proving convergence for the spiked TDP oscillations. 
What we will prove is that under the conditions stated above T(d\ ) can be written as 



h— +00 



l(d 1 )= ]T a 



(83) 



where 



a<*> = sl4>{s k )e s * d * / (dA/ds) \ s=s 



(84) 



The convergence of the sum in ( J83| ) is not uniform at d\ — > 0, but is uniform for any interval of d\ bounded away from 
0. We give the details of the relationship to the physical problem of Sec. Ill after we prove the main result above. 
The proof is organized with a set of lemmas. 

Lemma 1: Let s™ ld = s n + iir/2L and let H n be the horizontal path (as shown in Fig. |l4|) in the s plane from the 
imaginary s axis to the point s™ %d . Let 



'-H n 



F(s) ds 



be the integral of F(s) on this path, then N can be chosen so that for n > N 

\X H I < COnSt X n 2-p-(3d 2 /2L) _ 



(85) 



(86) 



Proof: The discussion of the roots of the spiked TDP Wronskian in Appendix I can be applied to the roots of A. 
For large n, the n dependence of the roots takes the form 



^(•Sra) = nir/L + const. + • • • 
K(s„) = -(3/2L) ln(ra) + const. + 



(87) 
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It follows that a constant can be chosen so that |s n | > const x n, and hence \4>(s)\ < 2irK/n p , for some constant K 
independent of n. Since 5ft(s) < on H n and d\ is nonnegative, we have that \e sdl \ < 1. With s = s n + iir/2L + u 
and with a running from — 3?(s„) to on the path H n , the integral must satisfy 



-«(..) df7 



where D(s) = e 2 A(s). On H n we have that 



and D(s) can be written as 



D(s) = P 1 (s n )e 



„2Ls 



-e 2is "e 2L<T 



where 



Pl(Sn) 



4 P 2(Sn) 



(89) 



(90) 



(91) 



(92) 



For n larger than some N we can make both IZi and IZ2 arbitrarily close to unity. For n > N it follows that the 
magnitude of the sum in square brackets in ( 91 ) must be larger than the second term, and hence 



\D(s)\ > \P 1 ( Sn )\e 2L °e~^ d -\ll 2 \ 
By choosing N large enough we can make |Pi(s n )| and [R-2\ larger than n-independent constants, so that 

\D{s)\ > const x e ^<r e -U( s )d 2 = congt x e 2L* e -n( Sn )d 2e -*d 2 

We have then that 



(93) 



(94) 



Pirn I < 



< 



K' 

nP~ 2 ' 

K' 
nP- 2 ' 



„3?(s„)d 2 



,K(s„)d 2 



1 



Md 2 -2L) da 



2L-d 2 



It follows from ( |S7j ) that we can choose a constant so that |e^( s ") d2 | < const. /rv- 3d2 ' 2L >, and therefore that 

const 



PrrJ < 



„p+(3d 2 /2L)-2 ' 



(95) 



(96) 



which was to be proven. A similar proof shows that the same result applies to \In_ n |. 

Lemma 2: On the arc T n from s™ ld to s™ d with |s| = |s™ ld |, the magnitude of A(s) is larger than some constant 
that is independent of n. 

Proof: We write 



A(s) = P 1 (s) [! + *(»] 



(97) 



where 



$( S )Ee 2s VP 2 ( S )/P!( S ) 



We must show that |1 + $| is bounded from below by an n— independent constant. We recall that s™ d = s„ + iir/2L 
and that $(s n ) = —1, so that 



-.mid 



P 2 «" d ) Pi(s„) 



^ 2 (s„) Pi(sr d ) ' 

For n larger than some N m - m we have that <j>(s™ ld ) is arbitrarily close to unity, so that 



(99) 
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*(C ) = 1 + A, (100) 

and Nmin can be chosen large enough to make |A| arbitrarily small. 
On the arc T n we write $(s) as 

<f>(s) = e 2sL s 3 (A 2 /A 1 )[l + p(s)} . (101) 

By choosing AVin sufficiently large, we can bound the magnitude of p to be less than an arbitrarily small n-independent 
constant. For s on the arc r„, we now write s as 

s = R n ie w , (102) 

where R n = |s™ d | and 9 is the counterclockwise angle from the positive imaginary s axis to the point s on r„. We 
denote by 9 n the angle to s™ d , so that s™ ld = R n ie • 71 and we write 9 = 9 n + 69. In terms of this notation we have 

*(«) ^ F F 1 (S9)F 2 (S9) , (103) 

where 

F ^(l + \)(l+p(s))/(l+p( S r d )) (104) 

— e 2iLR n [cos{0 n +86)-cos(9 n )] g 3i5e (105) 

F 2 = e 2L ««[ sin ( 9 «)- sin («"+' 5e )l . (106) 

The complex phase of $(s) is near zero at s = s™ ld , and decreases as s moves counterclockwise along the arc T n . We 
use the expressions above to find at what value 89* of 89 the phase of $ first becomes —ir/2. We note that \Fi\ = 1 
and on the top half of the arc F 2 < 1. The value of 89* is given by 

2LR n [cos((9 n + 89*)- cos(0 n )] + 389* + ({6 n + 89*) = -tt/2 , (107) 

where £(#) is the phase of Fq. We note that 

| cos(0 n + 89*) - cos(6»„)| > sin(6» n ) 89* . 



From (87) we know that sin# n decreases with n as ln(n)/n, and R n increases as n. Thus the first term in (107) is 
larger than the second by a factor that increases as ln(n). The third term, the £ term, decreases with increasing n, 
and we can bound it to be smaller than an arbitrarily small constant by choosing N min sufficiently larg e. Fr om these 
considerations it follows that we can choose iV m i n large enough that the magnitude of the first term in (107) is larger 
than, say, 2/3 of the magnitude of the left hand side, and hence 89* satisfies 

2LR n [cos((9„) - cos(6>„ + 89*)] > tt/3 . (108) 

Let us also take iV m i n large enough so that 9 n + 89* < ir /4. In that case we have 

sin(0„ + 59*) - sin( 



cos(6»„) - cos(6»„ + 89*) 



> 1 . (109) 



From this it follows that 2Li?„[sin(6' rl + 89*) — sin(6* I1 )] > 7r/3, and hence F 2 < e~ n / 3 . Since the deviation of |Fo| from 
unity is arbitrarily small, let us use |Fo| < e*'/ 12 and conclude that |<i>| < e _7r / 4 at the point along r n at which $ first 
becomes purely imaginary. As 89 increases, the magnitude of $ continues to decrease. It follows that for every point 
along the top of the arc 

|1 + *| > 1 -e^ 4 . (110) 



A similar analysis starting at s_„ shows that (|110j ) holds also for the bottom half of the arc. 
Lemma 3: On the arc T n from s™ id to s™ d with \s\ = \s™ id \, the integral 

1 f sms)e^+ d ^ 
IrM =2^J AW ( } 



satisfies, 
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/ , m const , . 

frMl < nP+i3d2/L) „ 2 , (112) 

where the constant is independent of n. 

Proof: On T n we have that 5R(s) < 5f(s„), so that 

| e sd2| < |gK( S „)d 2 ( 113 ) 



and, for some constant, the right hand side of (113) is less than const/n^ 3d2 / 2i - 1 . We have seen that |A| is bounded 
from below on r„. With the bound on \4>(s)\ from (jsil), we have then that 

Since the integrand is everywhere nonnegative, we have that 

e R to dl |tis| < / e u{s)dl \ds\ , (115) 

•/ arc 

where the arc is the half circle with \s\ = \s™ ld \ — R n in the left half plane. But, the integral along the half circle is 

/ e^'^ldsl =Rn / e- dlR '> sin8 d9 = 2R n e"*^ sine d(9 

J&rc JO JO 

< R n e- 2dlR " e ^d3 = £ (l - e -**») . (116) 
Jo «i 

This completes the proof of the lemma. 

Proof of main result: We define C ra as the integration path on the imaginary s axis from iQ(s- n ) — iir/2L to 
i$$(s n ) + m/2L, and we define 

l n {dx) = ( F(s)ds . (117) 

We let I n , c ioscd(di) be the integral on the closed path consisting of C n , of T n , of H n and of H_ n traced backwards. 
From the lemmas above we have 

\l n {dx) -Closed (di) | < |J r „(di)| + |X ff „(di)| + |Z ff _ B (di)|< ^(X/^)^ ■ ( 118 ) 
The integral on the closed path is 2m times the sum of the residues inside the path, 

n 

2n,closcd(dl) =J2 a(k)eSkdl ' ( 119 ) 
— n 

where a^> is the residue of s 2 (f>(s)e sd2 / A at s — s&. Since the only singularities of the integrand in the finite s plane 
are simple poles at the roots of A(s), these coefficients are those defined in (p4[). We have then 



— n 

and our main result follows from the fact that I{d\) is the limit of I n {d\) as n — > oo. 

We now turn to the role played by the Cauchy data. In the Green function solution for the waveform [see the 
discussion following (^lj)], a function of s appears representing the integral of the product of /_ (s, a;) and the combi- 
nation J(x,s) = —ipo(x) — sipo(x). In the case of the TDP or spiked TDP, f-(s,x) = e sx . Let us suppose that the 
support of the Cauchy data is confined to the region xi < x < x\ The Cauchy data then enters the s integral through 
the function 

J(s) = f 1 J(x, s)e sx dx . (121) 
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If the initial waveform tpo(x) satisfies 



dP+ 



dxP+ 



< bo 



then from integration by parts, we have 



J Xo 



e sx ip {x) dx 



= -u:;e sx £iMx)]dx 



SP+ 1 1x2 



x i „sx d p+1 
dxP+ 



T [tpo(x)] dx 



± 



sP +1 



For K(s) < the factor e s ( x X2 > in the last integral is < 1, so that 

\b Q {xx - x 2 )\ 



e 8X il) (x) dx 



< 



\s\p+ 1 



(122) 



(123) 

(124) 
(125) 



(126) 



If in addition to the constraint in (122) we have that the p th derivative of ipo{ x ) is bounded, then by a very similar 
argument we can show that </>(s), defined as e~ SX2 J'(s), satisfies 



* < 



const 

w 



(127) 



We can now apply the above mathematical results to the Green function integral from Sec. III. The waveform is 
given by the following integral along the imaginary s axis: 



ip(t,x) = 



1 

2i7T 



s(t—x) 



J(s) 



W{s) 



ds 



(128) 



where W(s) is given in (|57]) - ( |60| ) and has has the form W(s) = s 
and d8(i|). We can therefore rewrite the solution as 

ip(t, x 



-2 n -2sL 



A(s) in which A(s) is a special case of (T7 



1 



s 2 e s(t-x+x 2 +2L) 



A(,s) 



■ ds 



(129) 



The above proof requires that p, d%, d 2 and L satisfy the inequalities that follow (pT|) . The details of the proof show 
that the rate of convergence depends on these parameters. In particular, on the value of 7 = p+ (3c?2/2L) — 2. A small 
value of this parameter means slow convergence. We can now relate the details of the proof to the examples presented 
in Sec. Ill, and examine the interesting nature of the convergence of the seri es g iven there. We start by noting that 
a str aigh tforward computation of (f>(s) = e~ SX2 J(s), for the Cauchy data of jl9|), shows that p = 1. Comparing ( |78| ) 
and (|129() we see that 



d\ + d,2 = t — x + X2 + 2L 

= t — x + xi — x\ + X2 + 2L 



(130) 



An "obvious" choice is to take d\ = t — x + X\, the retarded time from the start of the reception of signals from the 
Cauchy data. At any x this is the equivalent of t — t m m- With this choice we are left with d 2 = x 2 — x\ + 2L. (Note 
that 2L — d 2 — x\ — x 2 is positive, as required in the proof.) The value of 7 for this choice is given by 



7 =P + 



M2 
2L 



2-3 



x\ - x 2 
2L 



(131) 



In our examples, we choose x x — x 2 , the range of support of the Cauchy data, to be 6 and we have L = 9, so 7 is 
unity. Suppose, though, that we had chosen x\ = 1 (as in our standard sine wave data of Sec. Ill) but had decreased 
the value of x 2 below our standard choice -5. The details of the proof show that convergence would require more 
terms for a given level of accuracy and that the series would fail to converge for x 2 < 11. This limit can be extended 
if we use initial data with one or more continuous derivatives, i.e. if p > 1. This rather unusual feature was, in fact, 
exactly what was observed in numerical experiments. 
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We point out next that Lemma 3 shows that convergence is not uniform in d±. As d% gets smaller, more terms in the 
series are needed. With our choice of di to be t — £ m ; n , this implies that the convergence of our QN series in (^3|) is not 
uniform as t — t mul — > 0, contradicting our claims of uniform convergence made following (^|). But note that we can 
choose d\=t — i m i n + 1, so that convergence is uniform for t > i m i n . In this case we have 7 = 2 — 3(xi — X2 + 1)/2L. 
For both our sine wave Cauchy data this 7 has a numerical value of 5/6, and the series is convergent. It is clear that 
that there is an interaction between the allowed range of the Cauchy data, and the range of t for which the QN scries 
converges. By moving the left edge of the support of the Cauchy data to the left by some amount S, we increase by 
35/ (2L) the value of t at which the series first converges. It should also be noted that the dependence on d\ explains 
why the QN series converge more quickly at early times than at late, a feature evident in Fig.KJ. 
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